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A particle subject to successive, random displacements is said to execute 
a random walk (in position or some other coordinate). The mathemati- 
cal properties of random walks have been very thoroughly investigated, 
and the model is used in many areas of science and engineering^^ as 
well as other fields such as finance 3 and the life sciences 4 . This letter de- 
scribes a phenomenon occurring in a natural extension of this model: we 
consider the motion of a large number of particles subject to successive 
random displacements which are correlated in space, but not in time. 
If these random displacements are smaller than their correlation length, 
the trajectories coalesce onto a decreasing number of trails. This sur- 
prising effect is explained and quantitative results are obtained. Various 
possible realisations are discussed, ranging from coalescence of the tracks 
of water droplets blown off a windshield to migration patterns of animals. 

The phenomenon we have discovered is illustrated in figure El It shows the 
positions x of 20 particles as a function of time t. Each particle undergoes a random 
walk, but the displacements of nearby particles are correlated, such that particles 
which are very close experience almost identical displacements. The paths of the 
particles are seen to coalesce onto a decreasing number of trails, along which the 
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particles follow almost the same trajectory. This effect is surprising because random 
walks normally reduce rather than accentuate inhomogeneities in the density of 
particles. In the following we explain this effect quantitatively and suggest some 
contexts in the physical sciences in which the effect can certainly be observed, and 
areas in biology where it may find interesting applications. 

We consider a system of particles labelled by an index i, having positions Xi(t) 
at time t. The particles are subjected to random displacements at times which are 
integer multiples of a small increment St. At time t = n8t the displacement is given 
by a random function f n (x) evaluated at the position Xi(t): 



The random displacements f n satisfy (f n (x)) = and (f n ( x )fn'( x ')) = 5 nn ic(x — x'). 
Here 5 nn > is unity if n — n' and zero otherwise, and we write (A) for the ensemble 
average of a quantity A (that is, the average over different trials of the random 
process generating A). The correlation function c(X) is even in X = x — x', decays 
rapidly as X — > oo, and is small when \X\ ^> £. A suitable choice (adopted in figure 
[Hand in the following) is c(X) = e 2 exp(— X 2 /2£ 2 ), e is the typical magnitude of each 
displacement and £ is the correlation length. Figure ^ used the values e 2 ~ 8 x 10 -6 
and f = lfr 1 . 

As an example of a situation where the effect could be observed, consider the 
motion of liquid droplets on a surface, moving in one direction under a constant 
force (rain blown off a perspex windshield is an example of this situation). If the 
surface is randomly contaminated, the wetting angle will be different on opposite 
sides of each drop, and the trajectory of the droplet will be randomly deflected. We 
are concerned with the case where the surface contaminants are smeared over an 
area large compared to the droplets (perhaps resulting from cleaning the windshield 
with a waxy polish), so that nearby droplets are deflected in the same direction. 
There is no interaction between the drops unless they are close enough to combine 
due to surface tension: we stress that the coalescence is that of the paths taken by 
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different drops, not of the droplets themselves. We model the motion of a droplet 
across the surface by a particle of mass m. At position r = (x, y) on the surface, 
the drop is subject to a force F(r) + F j, where F is the magnitude of a steady 
force acting in the direction of the unit vector j defining the y-axis and F(r) is a 
homogeneous and isotropic random force with correlation length £. We assume that 
the particles are subjected to a viscous resistive force proportional to their velocity 
across the surface, such that the equations of motion are 



(where p is the momentum of the drop). When the fluctuating force is weak, the 
trajectories are locally approximated by straight lines, with x approximately con- 
stant and with y increasing at a rate v y = Fq/tw). Furthermore, when the damping 



v y j — F(x(t),v y t)/mj. Under these circumstances the displacement x(t) perpendic- 
ular to the steady force satisfies 



When v y is large, the velocity of the particle fluctuates very rapidly, and in this limit 
it is reasonable to model equation (j3J) by the stochastic equation (JTJ). 

Why do the trajectories in figure ^ merge? The simplest argument is based 
upon linear stability analysis. The separation 5x(t) of two nearby trajectories varies 
exponentially in time, with rate A = t^ 1 (\og\8x(t)\/\5x(0)\) = (log |1 + f' n \)/5t . 
Provided e/£ is small, the magnitudes of the derivatives f' n = df n /dx are small 
compared to unity. Taylor expansion of the logarithm gives A ~ —^(f^ 2 )/5t, so 
that A is negative. In this case, most nearby trajectories approach each other with 
an exponentially decreasing separation, implying coalescence. This argument also 
indicates that if e/£ is large, the coalescence effect disappears: if the function /„ 
has a Gaussian distribution, the exponent A becomes positive when (f^ 2 ) exceeds 





7 is large- the velocity is given by the local value of the force: v(t) ~ p(t)/m 
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A more complete and concrete understanding of the coalescence effect is ob- 
tained by considering statistics of the density of particles, g(x, t) = J2i3(xi(t) — x). 
Translational invariance implies that an initially uniform density remains uniform, 
(g(x,t)) = Qq. Path coalescence is revealed by the density- density correlation 
function JC(x, x';t) = (g(x,t) g(x',t)) — Qo6(x — x'). Because of translational in- 
variance, the correlation function K, is a function of X = x — x' only: we write 
JC(x,x';t) = K(x — x',t). A tendency for particles to cluster is demonstrated by 
K(X,t) becoming large for X small, in the limit t — > oo. When the typical mag- 
nitude of the jumps is small compared to the correlation length £, we find that 
the correlation function satisfies a generalised diffusion equation, or Fokker-Planck 
equatiorP"^ 

The diffusion constant D(X) = [c(0) — c(X)]/St approaches zero quadratically at the 
origin: D(X) ~ kX 2 for X < £, where K = -\c"(fi)/5t. When X > £ the diffusion 
constant approaches a constant value, Dq = e 2 /5t. Equation (j3J) was obtained from 
the stochastic model, equation (JTJ), but in the appropriate limits the density-density 
correlation function of (j3J) also satisfies 1 ^ equation (j3J). 

Now consider the properties of solutions of equation (jlj). We note that equation 
(|3J) is in the form of a continuity equation, dK/dt + dJ/dx = 0, so that the integral 
of the correlation function over all X is a conserved quantity. The flux of the 
correlation function passing the separation parameter X at time t is 

J(X,t) = -^[D(X)K(X,t)} . (5) 

Consider an initially uniform distribution of density, with value £>o (corresponding 
to K(X, 0) = Qq). For X <ti £ the diffusion constant is an increasing function of X. 
Together with (jSJ) this implies an initial flux of correlation towards X = 0. At large 
times, K(X, t) is thus sharply peaked at the origin. For X ^> £, on the other hand, 
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the approximate solution of (JH) is 



K{X,t) ~ g 2 erf 



( 



X 



) 
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Using the fact that K(X, t) satisfies a conservation law, we deduce that the average 
number of particles condensing into each trail at time t is 



When t/5t is large, and X <C £ (but not too close to zero) the flux J(X, t) is found 
to be approximately uniform. This implies 



The 1/X divergence of (jHJ) is non-integrable, so that this expression must fail near 
the origin. An exact calculation shows that K(X,t) is finite at X = 0, the height 
of the peak is K(0,t) = exp(2nt). In summary, at large times the correlation 
function K(X,t) develops a correlation hole for ( « I < y/ADot reflecting the 
path coalescence effect. The results of numerical simulations shown in figure El 
confirm these predictions. 

In the remainder we discuss a number of possible realisations of the path 
coalescence effect. Unlike the case of liquid droplets described above, these are 
speculative, but they show that the path coalescence effect is likely to have a very 
broad range of applications. 

There are many potential applications in the physical sciences involving fast 
particles interacting with a random potential. Two examples of this type are the mo- 
tion of rocks rolling down a scree slope, and the motion of highly energetic electrons 
through disordered crystalline solids. There are also possible connections between 
the coalescence effect and a 'streaming' observed in the flow of electrons away from 
a constriction in a two-dimensional electron gas with very low scattering^, and we 
shall discuss this case in some detail. The experiment shows regions of markedly 
increased current density persisting to some distance from the constriction. This 
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was explained by showing similarity to simulations of independent electron motion 
in the smoothly varying random potential of the doping atoms. This model is es- 
sentially (J2J), with 7 = 0. Theoretical discussion of this systenP has emphasised 
that caustics are important in understanding the empirical results. We remark that 
these experiments might show an even more pronounced effect if dissipation were 
introduced. Equations (J2J) can show various types of behaviour; figure El shows nu- 
merically computed trajectories for one choice of parameters. These trajectories 
initially show evidence of the streaming effect discussed in references"^ with the 
associated caustics. For larger times the trajectories coalesce, in contrast to the 
dissipationless model where the streams eventually disperse. We note that dissipa- 
tion of the electron motion could be increased by increasing the temperature of the 
system. 

There are also potential applications in the biological sciences, involving the 
movement of organisms in response to small random fluctuations in their environ- 
ment. Our model provides a mechanism through which large numbers of organisms 
can congregate without communicating. 

One example of this type is the migration of animals across a nearly homoge- 
neous smooth terrain. Thus figure ^ could be thought of as a map showing paths of 
animals on an Eastward migration. The paths of the animals will be deflected by 
small random fluctuations of topography or vegetation. Our calculations show that 
the animals can be drawn together onto the same paths, even if there is no com- 
munication between them, and no gross features in the terrain favouring particular 
routes. 

A second example applies to simple organisms such as plankton which can 
move in response to changes in their environment, such as nutrient concentration. 
In cases where there are small, spatially correlated random fluctuations of the nu- 
trient concentration, the path coalescence effect could lead to unexpectedly large 
concentrations of organisms. Such a mechanism could be utilised by evolution, en- 
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abling simple organisms which cannot communicate directly to congregate for sexual 
reproduction. 
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Figure 1: Positions Xi(t) for 20 particles performing spatially correlated random 
walks. The trajectories coalesce. 



Figure 2: Statistics of the particle density g(x, t). a The density-density correlation 
function K(X,t) of the motion illustrated in figure ^ derived from simulations (o), 
compared with the limiting theoretical forms (JHJ), full line, and (jSjl. dashed line. 
Parameter values: e 2 « 1.25 x 10~ 8 , £ « 6.4 x 10~ 3 , and t = 5x 10 5 5i. The inset 
shows numerical results verifying that K(0,t) = exp(2ftt) for the same parameter 
values, b The mean number of particles in a cluster, N(t), for the process shown in 
figure ^ The results of the simulation (o) are compared with theory (JJJ), full line. 
Parameter values: e 2 ~ 2 x 1CT 7 , £ ~ 1.6 x 1CT 3 . Particles are considered to be part 
of a cluster of N if all of their positions are within an interval of length £. 



Figure 3: Trajectories (red) for damped motion in a smooth random potential 
described by equation (j2J. Particles are introduced with uniform density on the left; 
here the x axis is vertical. The force is a sum of two terms: a random force derived 
from the gradient of a potential, F(r) = — W(r) with (V(r)V(r')) = e 2 exp[(r — 
r') 2 /2£ 2 ], and a steady force F j acting in the y-direction. The potential is shown 
in green, higher values correspond to darker colours. Parameter values: F = 1.5, 
e 2 = 1.25 x 1(T 4 , f « 6 x 10" 2 , m = 1, and 7 = 4. 
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This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0303635v2 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0303635v2 
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